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This paper provides a preliminary assessment of electromagnetic wave 
propagation in focusing media which departs from those previously studied, 
ideal lenses and continuous media with square-law index variation. New 
approximate methods are described for obtaining the transverse beam width, 
phase constant, and ray trajectory in continuous lens-like media, and for 
determining stability conditions in lens waveguides, where the lenses con- 
tain splierical aberration and the continuous media contain fourth-order or 
higher-order terms of variation in index of refraction. 

It is proven that only in aberrationless-lens waveguides or in a con- 
tinuous medium with square-law index variation will the shape of a beam 
injected off axis or with an angle to the medium's axis remain constant 
about a beam axis which oscillates about the axis of the medium. In non- 
square law media the beam will spread, but knowledge of the coefficients 
describing the medium and the position and angle of the injected beam en- 
ables one to specify the maximum radius within which all of the energy 
will be confined. 

The following is an example of the type of solution obtained for non- 
square law media: for a medium characterized by the transverse index varia- 
tion 

n = n a (1 — §a 4 .r 4 ) 

ami assuming no index variation exists in the direction of propagation, the 
radius to the 1/e point infield is approximately 
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where the free-space wavelength X = n a \ and m is the order of the mode, 
m — 0, 1, 2 • • • . Quite generally, non-square law media show dispersion, 
and unlike the square-law media the various modes travel with different 
group velocities. Expressions are given to allow these effects to be evaluated 
for small perturbations on a square-law medium as well as for higher-order 
index variations. 

The transverse beam shape associated with any law of index variation 
is shown to be as well approximated (in the region of significant power 
density) by a cosine function or Gaussian function as is the field for an 
ideal lens-waveguide approximated by a Gaussian function in the presence 
of typical diffraction losses. 

Normal mode shapes are obtained for resonators with fourth- order and 
eighth- order mirrors by the method of Fox and Li; diffraction losses for a 
few Fresnel numbers are also given. In a certain range of Fresnel numbers, 
fourth- and eighth-order mirrors give lower diffraction losses than spherical 
mirrors. 

An approximate method for solving the paraxial ray equation for rather 
general (non-square-law) media is outlined. Requiring only reciprocity 
and symmetry about the medium's axis, it is shown that the radial position 
of the ray (x) is related to distance (z) along the axis of the medium by 

x = 26m cos m(3z 

where m = I, 3, 5,7 • • • . Moreover, it is shown that this series converges 
very rapidly, making it possible to get a good approximate representation 
with only a few terms. For example, for the fourth-order medium described 
above, an approximate solution is 

x = x {0.959 cos fa + 0.041 cos 3£z} 

where /3 = xq y/ 1.44a 4 . It is characteristic of all non-square law media 
to have a ray period 2ir//3 which is a function of the peak ray displacement, 
x . 

Lens waveguides with fourth- or higher- order terms in the focusing or 
index function can of course exhibit increasingly strong focusing for energy 
departing farther from the medium's axis, and if the lenses are suitably 
spaced might be useful in reducing the magnitude of beam wander due to 
imperfections or guide-axis curvature. However, for a given beam spot size, 
non-square-law lenses must be placed closer together than square-law 
lenses. 

Use of non-square law lenses or distributed media in transmission systems 
will most probably require repeater-system techniques which are operable 
with multi-mode signals at the receiver input. 
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I. INTRODUCTION 



Although a great deal of work has been done to describe electromag- 
netic wave guidance using a sequence of ideal aberrationless lenses, little 
has been done with more general lens-like media. The objective of this 
paper is to provide a beginning understanding of what happens to the 
important descriptive parameters of a light waveguide — wave phase 
constant, wave spot size, ray trajectory, stability restrictions — for 
general lens-like media, either continuous or formed from a sequence of 
focusing elements. Exact solutions for these parameters are difficult or 
impossible, a possible reason for little having been done.* A second result 
presented herein is a series of analytical techniques for obtaining useful 
approximate solutions representing a broad class of lens-like media; 
some of these techniques are most clearly presented by giving examples, 
which unfortunately leads to a rather large number of equations. To aid 
the reader in finding the section dealing with a particular topic, an out- 
line and brief resume is given below. 

As pointed out by J. W. Tukey, there is no a priori reason why ideal 
aberrationless lenses are best for use as a communication medium. We 
would like to know what does happen to wave guidance as we depart 
from aberrationless lenses, which is the case studied extensively in 
earlier work. Gas lenses have nearly constant focal length but their 
"principal planes" are actually curved surfaces. 1 E. A. Marcatili and 
D. H. Ring have pointed out that these curved principal planes have 
effects similar to those expected in plane lenses with spherical aberration. 
Thus the present work relates to existing gas lenses as well as to the 
question of whether or not to attempt to create a different form of 
inhomogeneous medium for light-wave guidance. 

Some of the present work dates back more than a year; impetus to 
putting it on paper was given by recent calculations planned by E. A. 
Marcatili and D. Marcuse. 2 These calculations showed that ray optics 
can accurately predict the loss even when half of the beam falls off the 
edge of the lens. Many of the conclusions drawn here are based on ray 
optics but since the controlling transverse variations take place over 
distances that are large compared to the wavelength, ray optics should 
give a correct conclusion. 

The subject is developed in this paper in the following manner. Sec- 
tion II covers ray optics for aberrationless lenses of any thickness, in- 



* As this paper goes into print, S. J. Buehsbaum points out the existence of an 
exact solution for a particular non-square-law /(j) originally derived with refer- 
ence to a quantum-mechanical problem. J. P. Gordon will report on this in a later 
publication. 
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eluding the case of a continuous medium. A ray-optic method for de- 
termining stability conditions applicable (using results of later sections) 
with generalized focusing elements is also covered here. Proof that the 
transverse field distribution of a normal mode beam injected off the 
medium's axis (or which goes there due to curvature, or displacement of 
the medium) is preserved only when the index of refraction decreases as 
the square of the distance off axis is presented in Section III. Section IV 
develops the proof that in a wide class of continuous focusing media, 
the ray paths are representable by a series of odd-harmonic cosine or 
sine terms. A definition of the focal length of a segment of an arbitrary 
medium is covered in Section V. Section VI discusses a new technique 
for obtaining approximately the wave phase constant and spot size for 
generalized lens-like media, stated generally and illustrated for fourth- 
order and eighth-order media. A characteristic ray angle and a charac- 
teristic ray period is defined for the generalized medium. Calculations 
for field distribution and diffraction losses in resonators and lens wave- 
guides using non-square-law elements are derived in Section VII. These 
calculations are the only ones available on fields and losses in non-square- 
law resonators. Section VIII discusses a solution for ray paths in pure 
fourth-order, or a mixture of second- and fourth-order, media; the ap- 
proach for extension to other media is indicated. Some further discus- 
sion and acknowledgements are presented in Section IX. 

II. TRANSMISSION IN AN IDEAL LENS-LIKE MEDIUM 

The medium referred to as ideal is one in which the index of refrac- 
tion has the form 3 (see Fig. la) 

n = n a (1 — £a 2 .c 2 ) (1) 

where 

n a = index of refraction on axis, x = 
x = transverse dimension 

a 2 = a constant. 
The index is independent of z, 

dn/dz = 0. (2) 

We can obtain a solution for the path of a light ray in such a medium 
using the well known paraxial ray equation 

6 2 .r _ 1 dn /.,•) 

dz 2 n dx 

Here, x is the position of the ray at some value of z. Equation (3) is 
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Fig. 1 — Ideal lens-like medium, (a) — two-dimensional continuous medium 
with index of refraction independent of z and varying with x according to (1), 
(b) — light raj' paths in medium a for parallel input rays, (c) — light ray paths 
in medium a with non-parallel input rays. 

valid over a broad range of conditions provided only that the ray path 

makes a small angle to the --axis. For a medium of the form of (1 ) 



and with the restriction 



we can write 



dn 

— = —n a aox 
ox 



haoA « 1 



(4) 



(5) 



dx 



1 dn 



dz 2 n dx 



= —doX. 



(0) 



The general form of (6) indicates an exponential solution for x as a 
function of 2, and it can be verified that 
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x = A cos \Za 2 z + B sin \/a 2 z 

is a solution, A and B being constants to be determined. 
Using the boundary conditions at z = 0, 

x = n = input ray position 

d.v 



dz 



= r, = input ray slope. 



We have 



x = r, cos Va 2 z + r t -\/\/a2 sin va-» z. 



(7) 

(8) 

(9) 

(10) 



This is very interesting in that it corresponds exactly to the ray be- 



havior in a sequence of equally spaced aberrationless lenses. ' It is 
important that all rays have the same oscillatory period regardless of 
input displacement or input slope, as sketched in Figs. 1(b) and 1(c). 
It is likewise significant that the effects of input ray slope and position 
are separable with respect to x, subsequent ray position. 

If a 2 is negative, the cosine and sine of (10) become cosh and sinh 
and a divergent medium results. 

A short segment of the distributed medium characterized by (1) 
may be assigned an equivalent focal length, even when it is not a "thin" 
or "weak" lens. With reference to Fig. 2, the most general case is one 
in which the index on axis, n a , is different from that of the surrounding 



LENS 




Fig. 2 — Diagram defining focal length and principal plane for an ideal dis- 
tributed lens. 
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medium w . The path of input ray with zero slope and displacement 
r, can be obtained from (10). 

The general equation for the ray slope, derived from (10), is 



— = — n y/a.2 sin y/a% z + n cos \/a 2 z. 
dz 

At the lens output the displacement is 

r = r t cos -\/a 2 1 
and the ray slope is, from (11), 

rj = —n y/ai sin y/a 2 1. 



(ID 



(12) 



(13) 



Because (5) holds, the refraction at the lens output surface is (see 
Fig. 2) 



Hence, 



d = 
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aa 
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wo cot -y/a 2 1 


a-< 


w a va 2 



(14) 



(15) 



The distance from the focal point to the principal plane where an equiv- 
alent thin lens may be placed (Fig. 2) is 



_ n w^ 

a-> ii a y/ai sin s/a 2 1 ' 



(16) 



This result was excerpted from the present work and given previously 6 
with the approximation appropriate to gas lenses Wo = n a . 

With the continuous focusing medium according to (1), there is no 
limit to the strength of the focusing effect; a stronger index gradient 
merely confines the normal mode energy more closely to the axis. How- 
ever, when a series of segments of the distributed medium are used as 
a waveguide, with gaps of a homogeneous medium in-between, there is 
a cut-off effect which occurs if the lenses are too strong in relation to the 
lens spacing. This is a well-known effect with thin lenses, 4 and the limit- 
ing conditions were derived from wave optics for the medium according 
to (1) by J. R. Pierce' and E. A. Marcatili. A ray-optic derivation is 
given here for the medium (1 ) to show the method which can be extended 
to arbitrary focusing media, including fourth-order or higher-order 
terms in x, which has not yet been handled by wave optics. 
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With reference to Fig. 3(a), we assume an input ray of zero slope 
and a displacement r, at the center of the first lens, where the longitudinal 
axis Z\ = 0. At Zi = t/2 and restricting our attention momentarily to 
the region < \Za 2 t < t/2, the output of lens # 1 is 



r i = r, cos (\/a 2 t/2) 

r i = — r, \Za 2 sin (y/a* t/2). 



(17) 
(18) 



Taking the approximation appropriate to gas lenses, no = n a , the ray 
output of lens 1 intersects the axis at a distance 6/2 from the end of 
lens #1 (Fig. 3a). 



6/2 = | roi/ro/ |. 



(19) 



Using our knowledge that (1 ) is symmetrical about re = and reciproc- 
ity holds, we can construct the remainder of Fig. 3(a) with a ray maxi- 
mum r, at the center of the second lens. The value of 6 corresponding 
to this condition is 



6 = -7= cot Wo* t/2) . 
Va 2 



(20) 




Fig. 3 — Ray path at one maximum lens spacing, (a) — weakest lens case, 
(b) — second stability band (stronger lenses than 3a). 
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It is clear that this is the largest permissible ray displacement at the 
center of lens # 2, since a larger value there would represent displace- 
ment growth without limit as2-» » . We can show that larger values 
of b violate this stability condition by writing the expression for the 
displacement at the center of the second lens. The input to the second 
lens is 

slope Vi-i = r,)/ (21) 

displacement r i2 = r ni + brai . (22) 

Then, from (10) and using a new origin of the longitudinal axis z 2 = 
at the input to the second lens, 

r i 

x 2 = r i2 (cos \/a 2 z 2 ) + ~7= sin \Za 2 z 2 . (23) 

Using (21) and (22) 

— = cos (\/a 2 z) (cos (\/a 2 t/2) — b y/a 2 sin (\/a 2 1/2)\ 

Ti ' ' - (24) 

— sin (\/«2 £2) sin (\/«2 1/2). 

Making the substitution 

2 
b = —r- cot (V«2 t/2) + 8 (25) 



and evaluating (24 ) at z 2 = t/2 we find 



.r 2 



Zi2 = t/2 



= -[1 + 8 sin (\/a2 t/2) cos (-y/oa </2)] . 



(26) 



In the region < \Za 2 1/2 < ir/2 both sin (\/a 2 1/2) and cos (\/a 2 t/2 
are positive; hence for 8 positive, the ray at the center of the second 
lens is displaced more than the input ray and a divergent path is followed 
in such a sequence of lenses. For 8 negative the propagation is stable. 
For ir < y/a 2 1/2 < 3ir/2 another passband occurs with the limiting 
value of b having a ray path as in Fig. 3(b). This leads to the stability 
conditions 



b ^ 77- cot (V<M/2) (27) 



applicable in the regions of lens length 



2026 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 

rnr < y/az t/2 < jit +- 



(28) 



where n = 0, 1, 2, 3, etc. 

For other permissible values of lens length we derive the stability con- 
dition with reference to Fig. 4. The input ray at z\ = 0, center of lens # 1, 
has zero displacement and a slope of r,'. Fig. 4(a) shows the ray path 
for the limiting value of b when ir/2 < y/ai t/2 < ir, constructed using 
symmetry and reciprocity as before. This yields 



2 .— 

b = — —y= tan (vao t/2) 



(29) 



This gives a positive value for b because the tangent is negative. 
Following the lines indicated above we can derive an expression for the 
ray path in lens # 2 



%1 = COS (Va2 Zi) 



02 



sin (\/a 2 t/2) + br/ cos (-\/a<i t/2) 



(30) 



+ —j= cos (\/«2 t/2) sin (-\/«2 z) ■ 

V«2 




vy jw \y x-i 



(b) 



Fig. 4 — Ray path at another maximum lens spacing, (a) — weak lens case 
(b) — second stability band (stronger lenses than 4a) . 
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We want the slope at the middle of this lens, which is obtained by putting 

z-> = t/2 in 

-f- 2 = — \/«" sin ( \/a-, z) { —7= sin ( -y/a 2 t/2) 
az» "V da 

(31) 

+ br/ cos (V'os f /2) > + r/ cos (\/a2 */2) cos (\/a2 2) • 

Letting 

2 _ 

b = - —7= tan (VM/2) + 5 (32) 

Vfl! 

we obtain 

1 dc, = x _ gv /- gin (y/^ t /2) cos (Voit/2). (33) 

ri dZ2 

Z2 = */2 

In the region 7r/2 < (\/a2 </2) < tt, the sine term is positive and the 
cosine term is negative; hence a positive 5 indicates instability, a slope 
at the center of the second lens greater than the slope at the center 
of the first lens and (as may be verified in (30)) a displacement at 
z-> = t/2 adding to the subsequent ray divergence. 

Higher order passbands occur, one sketched in Fig. 4(b). Thus, in 
the region 

(nv + t/2) < {Vat t/2) < (n + 1)tt 

n = 0, 1, 2, 3, 4, etc., the permitted range of lens spacing b is 

b £ - -7=tan(v / o 2 */2). (34) 

VO-. 

Equations (27) and (34) are identical to the relations obtained from 
wave theory. 7 

The above method for determining stability of lens waveguides can 
be applied for arbitrarily shaped lenses, as long as symmetry about the 
z axis and reciprocity exist. Numerical integration can be employed 
where the lens cannot easily be represented by functions in closed 
form. 1 

III. WAVE BEHAVIOR IN INHOMOGE.VEOUS MEDIA 

We will now consider media which are of the form 

n = n a [1 + f( x )] (35) 
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where 

fix) = —\a&~ — \a&? — \a& ••• (36) 

in which a 2 , a 4 , ae • • • are constants which may be positive or negative. 
Positive values represent convergent focusing and negative values 
represent divergent focusing. This represents a general medium re- 
stricted only to symmetry about the axis x = and uniform with 
respect to the direction of propagation. 

The normal mode for such a medium is characterized by the index 
variation. We will presently discuss the field shapes and losses for 
several special cases of (35); for now it is sufficient to note that a differ- 
ent equation for the index n means a different normal mode. 

An interesting general conclusion can be drawn by expanding the 
index n of (35) about some off-axis radius ri as in Fig. 5: 

n(x — ?'i) _ n(x ) _ nix = r t ) 
n a n a n a 

-f- x'{ —a 2 ri — 2a 4 ri' — 3a 6 ri ■ • • 1 



(*') 



2 -, r 4 



H <j— | —a-> — 6a 4 ri — 15a 6 ?'i - - , .,- 

( -'V 3 
+ ¥LL { - 12a 4 ri - 60a 6 ri 3 • • • I 
6 

+ ^- {-12a4 - 18(W---} 



in which 



X = (X — 7*1 ). 



With the approximation 



/(*) « 1 (38) 

the first term of (37) is unity. Then for a 4 , o 6 , and all higher order 
coefficients equal to zero, (37) becomes 

n(x) = n a {\ — aznx — ^ch(x) 2 }. (39) 

This is sketched in Fig. 5. Thus, for the medium according to (1), 
the effect of entering off axis at x = n instead of at x = is to intro- 
duce a term in the index which is linear in x. This term has the effect 
of delaying or advancing every region of the transverse cross section 
an amount proportional to the displacement from the axis; this term 
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Fig. 5 — Expansion of n(x) about point x = n . 



acts the same as a dielectric wedge and tilts the wave front. The path 
of any ray is described by ray optics in the continuous medium accord- 
ing to (10). The remaining term of (39) is identical to the original 
equation at n = 0; hence the normal mode for a wave entering at x= r\ 
is the same as at x = 0, with a direction change superimposed.* Thus, 
the beam follows a path sketched as in Fig. 6(a); a pure mode intro- 
duced into the square-law medium travels without change of shape 
if the change in position and direction of the beam axis is taken into 
account. If a plane wave is introduced, Fig. 6(b), the field is concen- 
trated periodically at the points where the beam axis crosses the axis 
of the medium. 

We can also see immediately that these properties are not charac- 
teristic of any medium in which a 4 or higher order terms arc present. 
In the expansion (37) there are, in addition to the term linear in x , 
additional terms in (.(•')" and (x'Y brought in by o 4 , a 6 , etc., which 
change the normal modes in a manner dependent on the particular 
value of i\ . When any one or more of a 4 , a B , etc., are non-zero, the 

* This was previously proven for a sequence of ideal lenses independently by 
H. E. Howe and J. P. Gordon. 
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Fig. 6 — Ray path representation of beam flow in the ideal square law-medium, 
(a) — normal mode off-axis input, (b) — mode mixture off-axis input. 

modes for a region off the axis of the medium are different from those 
on axis and the rather simple and attractive situation sketched in Fig. 6 
no longer exists. It also follows that ray paths also differ importantly 
from those sketched in Figs. 1 (b ) and 1 (c ) . 



IV. RAY PATHS IN NON-SQUARE-LAW MEDIA 

When f(x) of (35) and (36) contain non-zero a 4 or higher order 
terms, the general shape of f(x) can be something like that of Fig. 5, 
even though a 4 or some other terms are negative yielding a defocusing 
tendency. We shall be concerned now with arbitrary values of at , 
a 4 , a 6 , etc., within the limits on x such that 

<0 for,>0 



df{x) 
dx 



> for x < 



This assures that all rays will be bent toward the axis at any x, per (3). 
We can observe several additional features of the general case: 
(1.) As a consequence of (40), any ray will monotonically go through 
a maximum and return to the axis, as sketched in Fig. 7. Due to sym- 
metry about the z axis and reciprocity, the curve x = g(z) will have 
even symmetry about z 2 , zn , etc. and odd symmetry about z 3 , z B , etc. 
The period of g(z), (z 5 — Zi), will depend on the coefficients ch,ai, 
in (36). 
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Fig. 7 — Ray paths in a generalized lens-like medium. 
(2.) It will be instructive to note that 

fix) = - 



(41) 



with m — » °° is the special case of (36 ) which corresponds to a step change 
in index at I x I = a:* 



f(x) — > and n = n a for | a; | < \ a 

- f(x) » and (n - 1 ) « (n a - 1) f or | x | -> | a 



(42) 



We will see that this is a convenient bounding condition on non-ideal 
focusing media. In Fig. 7, the dotted line represents the ray path for 
such a medium. 

(3.) As already noted, the paraxial ray equation (3) is valid. Since 
f(x) in (36) is to be kept very small compared to unity, the solution 
for the position of the ray as a function of z, x = g(z), is related 



or 



dz 1 n dx dx 



fix) = f ^ [g(z)\ dx. 



(43) 



(44) 



* We will want to continue to use the paraxial ray equation and the associated 
condition, /(j) « 1. We shall use the region < | x \ < a but allow | x \ to approach 
a so closely that f(x) >£> 0. 
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(4.) Due to symmetry about the z axis, 

— [g(z)l = at values of z where x = g(z) = 0. (45) 

dz 2 y 



Obviously, 



s« f)1 -° 



(46) 



periodically at maxima for </(z), which will occur at values of z midway 
between the places where g{z) = 0, Fig. 7. 

(5.) It follows from the preceding notes that the function g(z) must 
have even symmetry about the values of z for which (46) holds, and odd 
symmetry about the values of z for which g(z) = 0, a consequence of 
reciprocity and the known symmetry about x = for/(x). Hence, the 
most general ray trajectory will be of the form 



= 2 b m cos m/3z 



(47) 



where m = 1, 3, 5, 7, • • • etc. We will proceed presently to find a par- 
ticular solution of the form (47). 

In all functions (35) where a 4 or some higher a n is present, the period 
of the ray trajectory depends on the peak amplitude. That is, if the 
medium characterized by (35) starts at z = (as in Fig. 8) and if a 
series of rays enter parallel to the z axis but at different values of x, 
the ray trajectories will have different periods. If a 2 , a A , etc. are all 
positive, the rays farther from the z axis will be bent more sharply and 
the picture qualitatively will be as sketched in Fig. 8. However, if a 2 




Fig. 8 — Raj' paths in a generalized focusing medium for various parallel 
input rays. 
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is positive and a\ negative (other a n = 0) the rays entering farther 
from the z axis will have a period longer than those entering near the 
z axis. As previously noted, all rays have the same period (Fig. 1) 
only when all a n except a* are zero. 

V. FOCAL LENGTH OF A SEGMENT OF AN ARBITRARY MEDIUM 

Suppose a ray at radius x passes through a section of arbitrary medium 
(Fig. la) of length t. We can ascribe a focal length to each x position 
as follows: the phase difference Ap between a ray on-axis and a ray at x 
is given by 

&<p = _I n a t{%a 2 x 2 -}- §a 4 x 4 - • • }• (48) 



We require t in (48) to be so small that 

A^«^' (49) 

A 

and we represent the segment of focusing medium as a thin lens; for 
any thin lens the focal length / is 

(50) 



/ = 


irx 
AoA<p 


Then, combining (50) and (48), 




f — 


1 



(51) 
n a t { ai + (I4X- + deX* • • • \ 

If a? = 0, the focal length approaches infinity as | x \ — > 0. 

If the medium is represented mainly by the square law term a 2 with 
a small a 4 perturbation 

f = nja 2 (l + R) (52) 

where 



R = a -^. (53) 

The ratio (53) appears repeatedly in describing slightly non-ideal media. 

VI. NORMAL MODE PROPERTIES IN ARBITRARY FOCUSING MEDIA 

It would be nice to have solutions to Maxwell's equations for a medium 
according to (35) and (36) but this has not yet been achieved. E. A. 
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Marcatili has solved Maxwell's equations and found normal mode field 
properties for the ideal square-law medium (only a 2 present), 7 and for 
a perturbation thereon (small 04 in a medium principally characterized 
bya 2 ). 8 

In this section, an approximate method for obtaining some of the 
properties of the modes of an arbitrary medium (pure fourth order, for 
example) is described. This is of interest as guidance on whether or 
not to put the effort into getting better solutions. 

We use as limiting cases the solutions found by Marcatili when only 
a<i is present, and the known solutions for a step transition in dielectric 
constant, which as previously noted in connection with (42 ) corresponds 
to a very high exponent on the x n term of (36). These two cases can be 
viewed as limits to the possible solutions for arbitrary (36), within the 
bounds of (40), and useful results inferred about the intervening cases. 

We note that the transverse field distribution for the first-order mode 
with the step-change of index, Curve I of Fig. 9, is of the form 



E = cos — 



© 



(54) 



1 

n 

m 

Bf 






-V 



^ 






n = n a [i-(i/a)°°] 
n = n a [1 - 5.28 x 10 20 1 8 ] 
n = n a [1- 1.975 xio 7 i 4 ] 
n = n a [i-2i2] 
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Fig. 9 — Index of refraction versus normalized transverse position (x/a) for 
several media. 
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where the index change occurs at x = ±a. From Marcatili's wave solu- 
tion, 7 the transverse field shape for the lowest order mode in a medium 
characterized by 

n = n a (1 — 2 a 2# ) 



is given by 



where 



E = exp 



[- 



(55) 



(56) 



2 ~ Vt (o,)» 



(57) 



X = \o/n a 

\ = free space wavelength. 

We now match the fields given by (54) and (56) by setting (54) equal 
to (56) at x = x; this yields 

a e = 1.315 x (58) 

where a e is an "equivalent" or effective half -width for the square-law 
medium. The corresponding field shapes and plots of index variation 
are given by Curves I and IV in Figs. 9 and 10. 




Fig. 10 — Normal mode amplitude vs normalized transverse dimension 
(x/a) for media of Fig. 9. 
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As already indicated, these two cases are limits on the index varia- 
tion. It is inferred that an arbitrary index distribution (35) will have a 
field distribution similar to that of Curves I and IV of Fig. 10 when x 
is related to the effective width a e to yield an index variation approxi- 
mating Curves I and IV of Fig. 9. This idea has been carried through 
to specify a procedure for defining a e for arbitrary media, and to yield 
particular solutions for several specific index variations. Approximate 
phase constants for the various modes of arbitrary media are also 
found. 

To document the inference on normal mode field shape for the non- 
ideal medium, two specific cases were considered in a little more detail. 
These were represented by Curves II and III of Figs. 9 and 10. Pure 
fourth order and pure eighth-order index variations were taken, the 
constant a n in (36) being selected to give an "eye-ball" fit to Curves I 
and IV of Fig. 9; the actual numbers for (a„/2) are on Fig. 9. Then, 
using the computer-simulation of a resonator as first employed by Fox 
and Li, 9 the normal mode field distributions were computed. More will 
be said about these calculations presently, but for now we note the 
conclusion on the normal mode field distributions for the eighth-order 
and fourth-order index variations, Curves II and III of Fig. 10. The 
choice of a„ for these curves was a first guess, not an optimized choice. 
Nevertheless, the fields do correspond very well to the limiting cases, 
Curves I and IV. We shall see that differences between Curves I, II, 
III, and IV in the region x < 0.9a are comparable to the differences 
between the true fields in the system of ideal lenses and the Gaussian 
approximations commonly used to represent them. 

The effective medium width a e can be defined for a more general 
non-ideal medium as follows: the normal mode in the medium charac- 
terized by (35), (36) and (40) will have a shape similar to that of the 
associated square-law medium or step-change medium when a c is defined 
by 

// > (1.315) 4 ; <AY ,- Q , 

f(x = a e ) = -r-^-b^-J. (59) 

This comes from making f(x) of (35) for the arbitrary medium equal 
to r times the same quantity for an equivalent square-law medium, both 
at x = a r , and then using (57) and (58) to eliminate the a 2 of the equiv- 
alent square-law medium. The constant b is unity for the lowest order 
mode, and is shown in Appendix I to be 

b = li/^p- (60) 
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for the TEM m mode.* The constant r varies slightly depending on the 
exact form of f(.v); in the examples plotted in Fig. 9, r is 1.72 for the 
fourth order case and 1.43 for the eighth order case. When considering 
a medium which is very nearly square law, r may be taken as unity, 
and for fourth or higher order media a value r = 1.5 is a good value. 

Following known theory for metallic waveguides, the normal mode 
for the step-change index can be represented by two plane waves travel- 
ling at an angle a to the axis of the medium, where 

sin a = — (61) 

and X c is the cutoff wavelength in the guide containing everywhere the 
index n a . Restricting our interest to the region X/X c <K 1, the phase 
constant is given by 



fi = f Vl - sin* a = y U " W-'i- (62) 



Xc = r _ifL TT ( 63 ) 



— V 1 — sin- a = - 

A , 

Higher modes are accounted for by noting 

4q 

(m + 1) 

giving 

X { 32 \a/ J ( 64 ) 

X = \ /n a . 

Since the principal part of the field distribution has been shown to be 
essentially the same for all/(.r) restricted by (40) provided a e is defined 
as in (59), we can expect that the expression (64) will give a good ap- 
proximation for the phase constant in a medium characterized by any 
/(•r). 

We proceed to write down the expressions which result. 

For the medium described by 

n = n„ (1 - \a*£ - \a^) (65) 

the application of (59) leads to 

- - L315,, tuW \/\ {m 

* Following Fox and Li and Boyd and Gordon 10 the mth order mode has (in + 1) 
field maxima in the transverse cross section. 
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with R = aide/di . The ratio R is the ratio of the fourth order term 0.5 
atx* to the second order term 0.5 a 2 x 2 at x = a e , the equivalent half- 
width of the medium. Thus the expression (66) gives a e implicitly. 

When R « 1, we can rewrite (1 + R)' = 1 + $ B; under this con- 
dition, putting (66) into (64) gives for the medium characterized 
principally by a square-law distribution with small fourth-order varia- 
tion, (taking r = 1), 



2ir .- \ (m + l) 2 , 0.275(m + l) 2 a 4 X l 



(67) 



We can compare this result to Marcatih's direct wave solution 7 in the 
limit of a 4 = 0. The functional dependence on X and a 2 is identical; 
the coefficient of y/ai in (67) compares to the correct one as follows: 





Coefficient o 


Voi 


III 


Marcatili 7 


Eq. (67) 





0.5 


0.357 


1 


1.5 


1.02 


2 


2.5 


1.78 


3 


3.5 


2.6 


4 


4.5 


3.43 


m —* oo 


m 


0.892 m 



Thus, the approximate value of the phase constant is correctly given for 
all modes at any wavelength by the approximate theory outlined above, 
but the unique integral relationship for the phase differences between 
the modes in the square-law medium is not given. 

A comparison between (67) and some unpublished results Marcatili 
obtained using a direct perturbation solution of Maxwell's equations 
shows identical dependence upon a 4 , a 2 and X with comparable but some- 
what different constants. 

The above comparisons lead one to have confidence in other results 
obtained from (59) and (64) for which there is no previous information. 
For example, we can get the phase constant corresponding to medium 
(65) with a 2 = 0. Equation (66) leads to 



r I/6 (1.315) 2/3 ,2/ 3 X 1/3 
a < = ^7i b ^7-6 



04 

/-m -I- 9 K 

= 0.876 



( m + 2.5 \ 
V 2.5 ) 



1/3 .1/3 

04 175 



(68) 



for the equivalent half-width of the medium, and 
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for the phase constant. 
Likewise, for a medium 



( m + 2.5 V 
\ 2.5 ) 



0-f-0**(a*) w , m + s fr m 



n = n a (l - \a&*) (70) 

the application of (59) and (G4) gives 



and 






2* p.. 0.2 (m + 1) 2 

= 0.323X a s 



(71) 



Y~— - / m + 2.5 Y 4 - (72) 

For the limiting case of a step change in index 



,*-'\ (m + l y. < 73 > 

X 16 a- 

Looking at the change in /3 as the exponent in variation of index of 
refraction changes from 2 to 4 to 8 to <x>, (see (67) with o 4 = 0, (69), 
(72) and (73)), we observe that the /3 dependence on X goes smoothly 
from X° to X*, X ' 6 , and X 10 . It seems certain that the square-law medium 
is unique in having beats between modes being independent of X. 

There is a physical explanation for the increasing dependence on 
X for increasing exponent in the index variation. For the step change in 
refractive index the modes are all contained in the same transverse 
space, illustrated by Curves I of Figs. 11 and 12; by definition no energy 
can exist at x > a. However, with the square-law variation in refrac- 
tive index, Curve IV of Fig. 9, energy can and does exist at x > a 9 
and the penetration of the field at x > a e increases as the mode index 
m increases. This is illustrated for m = and m = 2 as Curves II in 
Figs. 11 and 12 respectively. Since the higher order modes occupy 
more transverse width for increasing m, X c decreases more slowly for 
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Fig. 11 — Normal mode amplitude versus normalized transverse dimension 
(j/o) for step change refractive index and square-law refractive index. 
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Fig. 12 — Normal mode amplitude versus normalized transverse dimension 
(x/d) for a higher-order mode in the step change index and square-law index media. 
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the square-law medium than for the step-change refractive index. That 
is the basic reason for the factor 6, (60) appearing in (59). Similarly, 
the increased field penetration at x > a e for the square law compared 
to the step-change index makes the angle a smaller for the square-law 
medium than in the step-change medium, and consequently the (3 is 
less altered from the value 2ir/X. As the exponent in the terms of f(x) 
increases from 2 toward °o the behavior in /3 and X c approach those for 
the step-change medium because the more rapid variation in index at 
x > a e prevents the field from penetrating that region as much. 

One can write explicitly an expression for the half -width w e of the 
lowest order mode's field for the various non-square-law media using the 
relation 



w, = 



1.315 



(74) 



where a e is again defined by (59). As a consequence of the method of 
defining a e , (74) gives the exact spot size for the pure square-law me- 
dium. For the medium which is principally square law with, a small 
fourth-order term, (65) with R < 1, 

w -=VlM i -°-^- (75) 

For the pure fourth-order medium, 



02" 



,1/3 

w. = 0.666 ^— (76) 



(i\ 



1/0 



and for the pure eighth-order medium, 



X - 2 
w e = 0.703 ^— ; . (77) 



0.1 



It is possible now to specify unique rays to consider characteristic 
of the particular mode in each of the media. These rays are the normals 
to the two plane waves which in combination give approximately the 
transverse field distributions in the manner outlined above. Each of these 
normals makes an angle a to the axis of the medium as in equation 
(61). Again confining our interest to waves far from "cut-off" 

X(m + 1) 



" ~ 4a e 
where a e is again defined by (59). For the medium (65) with R < 1 



(78) 
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and 



a ~ 1315 J 
ae ~ at* V 



m 4- 2.5 
2.5 



- { 1 - 0.55 

7T 



(m+ 1) 



ra-J-2.5~| Xa 4 i 
L^5~~ W»f """ 



2.97 



, /m + 2.5 L n ,, Xa 4 (m + 2.5\\ ' (80) 

This expression also gives the characteristic ray angles for the square 
law medium by letting a 4 = 0. 

For the fourth-order medium, (68) gives a e and the characteristic 
ray angle is 



a = 0.285A" W F 



,,, ,..„ (m + 1) 



I' m + 2.5 T" * 
. 2.5 J 



(81) 



For the eighth-order medium, (71) gives a e and the characteristic ray 
angle is 



n s n i ("* + 1) 

« = o.27ix°V r^ 



W 



+ 2.5T" 
2.5 J 



(82) 



Having a characteristic half -width a e and ray angle a for each me- 
dium, we can easily calculate a characteristic ray period for each me- 
dium. With reference to Fig. 13, the ray at an angle a with the z axis 
has a period X p 



, 4a e 

A„ — . 



;83 



For the pure square-law medium, using (66) and (80) 
4 X 1.315 X 2.97 



X„ = 



1 [ m + 2.5 1 _J 




a P i 




Fig. 13 — Diagram for determining an equivalent ray period in non-square- 
law media. 



LIGHT PROPAGATION 2043 

or 

. 8.8 (m + 2.5) (su x 

a 2 J 2.o(ra + 1) 

We can compare this ray period to that derived from the paraxial ray 
equation; that solution is (10) from which the ray period is seen to be 
2w/-\/a-2 . Equation (84a) corresponds quite well to 2ir/\/a-i for the 
lowest order mode (m = 0) including independence of X! Equation 
(84a) is also relatively independent of mode order, m. 
For the pure fourth-order medium the ray period from (83 ) is 

. 12.3 [ m + 2.5 f 1 ,om 

Xp " V%77-3 [-^5— J (m-r-1)" {bo) 

There is no previous theory for comparison and, as already noted, 
ray theory alone does not define a unique period; unlike the square-law 
case, it depends on the initial ray slope. 

VII. FIELD DISTRIBUTION AND DIFFRACTION LOSSES IN RESONATORS AND 
LENS WAVEGUIDES USING NON-SQUARE-LAW ELEMENTS 

In this section there are recorded new computations of normal mode 
field distributions and diffraction losses for a few resonators using non- 
spherical mirrors. The results are applicable to lens waveguides using 
non -square-law focusing elements. No other information is known to be 
available on these configurations. 

The motivation for the work was to determine the normal-mode field 
distributions for non-square-law continuous waveguides to support the 
approximations made in the preceding section. Thus, most of the con- 
ditions selected represent relatively close spacing of weak lenses. A few 
instances are cited where cutoff effects associated with too great spacing 
of lenses were observed. 

The resonator whose descriptive dimensions are given in Fig. 14 is to 
be used to represent the lens guidance system of Fig. 15. The equiva- 
lence is exact if one includes the absorbing screens in Fig. 15; if the dif- 
fraction losses in the resonator are not large very little effect is expected 
from omitting the absorbing screens. Fig. 15, in turn, can represent a 
continuous guidance medium when the focal length of the lenses is ap- 
preciably greater than the lens spacing. Then the beam size is veiy 
nearly the same at the lens and at the plane midway between lenses, 
and the electrically thin lenses of Fig. 15 can be representations of 
segments of the continuous medium of length s. For the square-law 
medium, this approximation is that the focal length expression (16) is 
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Fig. 14 — Definition of parameters in a resonator for use in Fox-Li type calcu- 
lations. 



well represented by 



f = 



Wo 

n a aii 



(86) 



when yfai t « 1. 

In representing a continuous medium by a resonator it is important 
to have a measure of the accuracy of the approximation. We obtain 
this and observe an interesting outcome by comparing the phase shifts 
along equi-power contours (i.e., contours within which the fraction of 
the total power in the beam is constant) for the resonator with spherical 
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Fig. 15 — Transmission medium analogous to Fig. 14. 



LIGHT PROPAGATION 2045 

mirrors and for the square-law continuous medium, 

n = n a (1 — \aix ). (87) 

For the square-law medium of length s, the phase difference for a ray 
following the axis and one at radius .r is 



A<p c = — s( -%a 2 x 2 ) . (88) 

Ao 

For the resonator with spherical mirrors the phase difference between a 
ray on axis and one following an equi-power contour is 



2 
By making 



A<* m ^(-2A). (89) 

Ao 



A=^-V (90) 

4 

it is apparent that (88) and (89) are identical when x = x . It may be 
shown that the beam spot size at the center of the resonator is the same 
as the spot size of the continuous medium when 

Wa = JhL JM (9i) 

where R c is the field curvature at the confocal spacing of mirrors asso- 
ciated with the mid-plane spot size w .* The equi-power contours are 
given by 

■L = ^ = vTT> (92) 

.To Wo 

« - & (93) 

w, = spot size at the mirrors. 
Using (91), (92), and (93), the value of both (88) and (89) is 

A^, = A^ = ^.r 2 * (94) 

Ao Kc 

Thus by choosing the continuous medium and spherical mirror system 
to have the same mid-plane spot size, the phase shift along any contour 
enclosing equal powers is identical, regardless of lens spacing.! 

* See Ref. 10. 

t Within, of course, the region of stability for the resonator. 
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For reference it may be noted that the mirror radius of curvature 
b and focal length / are related to R c and s by 1 

h > = 2f= (1 + J )Rc (95) 

2 

R c = —7= = V2sb' - s 2 . (96) 

Vfl2 

Useful forms giving the mid-plane spot size w and spot size at the 
lens w s for the lens system of Fig. 15 are: 

X (At* J\k 



W> = /|/«f (4/s-s)* (97) 



-i/l^i'-if (98) 



Wo 



w - = 4/ i t: a 5, 09) 



X (/«)' 



It seems certain from (94) that the normal mode field distribution 
for a continuous medium will be well represented by the normal mode 
of a resonator having the proper correspondence between A (Fig. 14) 
and the coefficients of (36), provided that the mid-plane of the resona- 
tor field does not differ appreciably from the field at the mirror. 

The proper value for A of Fig. 14, to represent a general medium (35) 
and (36), is obtained by taking x = x (Fig. 14) and letting A<p m = A(p e 
as was done in connection with (88) and (89). This yields 



s 



/l 2.1 4,1 6 \ 

t 5 aiX + 9 a * x + o " ' ) ' 



A = - ( - aw 4- - a\x 4- - a& • • • ) . (100) 

Selection of a mirror spacing suitably small can be done for the square- 
law medium with the aid of (92). For example, with 

ai = 4 (meters) -2 

s = 0.2 meters 

wo = 0.319 X 10 -3 meters 

R c = I meter 



one computes w 3 /w = VToi. This is based on the Gaussian function 
approximations for the actual fields. A computer determination of the 
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actual field distributions* has been carried out using the method of Fox 
and Li and the results are plotted in Fig. 16 for a Fresnel number N = 
1.38. t These are true normal mode field distributions for the resonator; 
some perturbation from the Gaussian shape is caused by the 5.27 per 
cent power diffraction loss on each reflection but this does not account 
for the fact that the midplane 1/e width is approximately 7.5 per cent 
less than the mirror field instead of 2 per cent less as predicted by (92). 
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Fig. 16 — Computed normal mode fields at the mirrors and at midplane for a 
spherical mirror resonator. 

See Fig. 17, where the diffraction loss is increased by reducing the mirror 
size; the mirror spot size changes very little. Also see Fig. 18, where the 
fields, both the real fields developed on the computer simulation and the 
Gaussian theoretical approximation thereto, are plotted. We see here, 
with N = 1 where diffraction losses are less than 0.1 per cent per re- 
flection, that the Gaussian approximation is too narrow at the mirrors 

* In all of the computer determinations of field distributions and losses, 100 
radial intervals were employed in representing the transverse field distribution. 

f In the terminology of this paper, N = r 2 /s\ where 2r is the mirror diameter 
and s is the mirror spacing. 
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Fig. 17 — Computed normal mode fields at the mirrors for the resonator of 
Fig. 16 with two mirror diameters. 

and too broad at the midplanc; whereas w s /wq = \/2 from the Gaussian 
approximation, the actual ratio of widths at 1/e is about 1.58. The 
prolate spheroidal wave functions give computed results in much better 
agreement with the Fox-Li machine computation ' 1 and we presume 
the latter to be more accurate than the Gaussian approximation. 

The transverse fields corresponding to the fourth -order index variation, 
Curve III of Fig. 9, is shown in Figs. 19 and 20. The mirror spacing, 
losses and Fresnel numbers arc listed on the figures. The midplane field 
corresponds as well to the mirror field as did that for the square-law 
case, Fig. 16, and the spacing s is judged adequately short. For the same 
Fresnel number, N = 1.38 and the same spot size, the fourth-order 
medium had 3.91 per cent loss per reflection compared to 5.27 per cent 
for the square-law medium. This is because the field at x = a e falls off 
more sharply as shown on Fig. 10. The field plotted in Fig. 10 as Curve 
III was obtained from the computer data which produced Fig. 20, 
with the reflection number selected to minimize the amplitude of the 
higher modes still present in the data for Curve I of Fig. 20. 
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Fig. 18 — Comparison of normal mode fields in a confocal square law resonator 
as determined by (1.) Fox -Li type calculations and (2.) Gaussian analytical 
approximation. 
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Fig. 19 — Computed normal mode fields at the mirrors and at midplane for a 
fourth -order mirror resonator. 
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Fig. 20 — Computed normal mode fields at the mirrors for the resonator of 
Fig. 19 with two mirror diameters. 

Figs. 21 and 22 show computed field distributions for a resonator 
with eighth-order mirrors chosen to match the medium of Curve II, 
Fig. 9, at the same mirror spacing used for Figs. 16-20, 0.2 meters. 
The value of A is computed from (100). In Fig. 21 we note the midplane 
field is quite unlike the mirror field, and in Fig. 22 we note the loss 
does not decrease smoothly for increasing mirror size, but remains more 
or less independent of mirror size. This is due to the fact that the focal 
length, which is a function of radial position on the mirror according to 
(51), is too small compared to the mirror spacing at the outer edge of the 
mirror. For example, for the Curve III of Fig. 22, the focal length is 
0.077 meters at the edge of the mirror. For larger mirrors than those 
represented in Fig. 22, it was found that both the loss and field distri- 
bution became constant, independent of mirror size. The fields striking 
the outer edges of these mirrors are radiated because they are reflected 
through the axis of the resonator at such a sharp angle as to miss the 
opposite reflector. Fig. 23 shows a rough sketch of this situation. (Note 
that the midplane fields resulting would be composed of a mixture of a 
propagating field and a radiated field.) An analogy can be made to the 
unstable region which occurs for all rays in a resonator with spherical 
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Fig. 21 — Computed normal mode fields at the mirrors and at midplane for an 
eighth-order mirror resonator. 




Fig. 22 — Computed normal mode fields at the mirrors for the resonator of 
Fig. 21 with three mirror diameters. 
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Fig. 23 — Diagram illustrating over-focusing losses in a resonator. 

mirrors when / < s/4 or / < 0.05 meter when s = 0.2 meters as in Figs. 
21 and 22. 

The corresponding situation in a lens guidance system is shown in 
Fig. 24. The energy near the beam axis is weakly focused leaving a rather 
large spot size. The energy extending to larger and larger radii is even- 
tually focused too strongly and misses the next focusing element, illus- 
trated by the dotted lines of Fig. 24. Thus, ordinary field spreading, 
causing the usual diffraction losses, is mixed with over-focusing losses 
in a certain region of the parameters for a non-square-law resonator. 
It is believed that this accounts for the loss differences tabulated on 
Fig. 22. 

The objective of determining the normal mode field for Curve II of 
Fig. 9 can be met by making the mirror spacing in the equivalent 
resonator smaller. This has been done in Figs. 25 and 26. In Fig. 25, 
the midplane field and mirror field are shown to agree quite well. In 
Fig. 26, some variation in field shape results with reduction in mirror 
size, but this is in part due to higher order modes still present in the 
field plotted for Curve I where the loss is extremely low. For the normal 
mode shape plotted as Curve II in Fig. 10, a careful selection of the 
resonator reflection number was made to minimize the higher-order 
mode amplitudes. 

Using very closely spaced mirrors with so little curvature over the 
area where most of the beam energy is located, as in Curve II of Fig. 26 




Fig. 24 — Diagram illustrating over-focusing losses in a lens guidance system. 
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Fig. 25 — Computed normal mode fields at the mirrors and at midplane for an 
eighth-order mirror resonator representing a shorter segment of the same medium 
as in Fig. 21. 
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Fig. 26 — Computed normal mode fields at the mirrors for the resonator of 
Fig. 25 with two mirror diameters. 
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Fig. 27 — Computed normal mode fields at the mirrors and at midplane for a 
plane parallel mirror resonator analogous to the resonator of Fig. 25. 

and the associate medium curvature curve II of Fig. 9, one might wonder 
whether the guidance is caused by diffraction as with plane parallel 
mirrors or really by the curvature of the mirrors. Fig. 27 shows the fields 
and loss for the plane-parallel mirrors corresponding exactly to Curve II 
of Fig. 26 for the eighth-order mirrors. The loss is <0.45 per cent for 
the eighth-order mirrors and 1.4 per cent for the plane parallel case. 
Evidently the guidance is importantly determined by the curvature. 



VIII. RAY PATHS IN NON-SQUARE LAW MEDIA 

In this section, the general properties of media described by (35), 
(36) and (40) are exploited to obtain approximate solutions for the 
paraxial ray equation (3) in several non-square-law media, and to indi- 
cate an approach which reduces to straightforward algebra the problems 
of getting similar solutions for other media. 

It has been shown that (47) gives the general form of the ray path 
for the media of interest, and it was noted in connection with Fig. 7 
that the limiting case of a step-shift in index of refraction results in a 
triangular waveform for the ray path. For a triangular waveform the 
coefficients b m in (47) are in the ratio 1, 1/9, 1/25, etc. for m = 1, 3, 5, 
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etc. For the square law medium, a simple cos /3z represents the ray path, 
as given in (7). Hence, we know that for any medium of interest the 
series (47) will converge very rapidly. We can hope to get a good ap- 
proximation with only a few terms of (47). We now outline the results 
of such a solution for a medium described by (65). We set 



x = Ci cos (Jz + c 2 cos 3/52 
which is equivalent to 

X = x |(\ + -J 

x = Xo{l + A sin &z\ cos /3z 



(101) 



A 

3P ) cos 0z — — cos 3/3z 



(102) 



with 



A = - 






h c i) 






A 



4 + A 



Ci -\- c 2 = x . 



(103) 



(104) 
(105) 



In form (102) it is seen that the maximum x is xo , occurring at z = 0, 
for any value of A ; this is a useful form in visualizing the effect of at as 
a perturbation on a medium mainly controlled by a 2 . 
Using (101 ) we find 



j2 

_ _ (x) = — Ci/T cos /3z — 9c 2 /3" cos 3/3z. 
dz- 



Also, from (65) 



-T- U(%)] = —aw — 2a>i% • 
ax 



(106) 



(107) 



We can get the solution to the paraxial ray equation in the form (43) 
by equating (106) to (107) with x replaced by (101). Equating the 
coefficients of the cos fiz terms yields 



= y/ a2 {1 + #(1.5 + 0.376,4 + 0.1875A 2 )} 5 



(108) 



where 



R = 



QiXo 

a 2 



(109) 
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We note that the restriction (40) requires (107) to be greater than 
zero, or 

R ^ ~h 

This is a limitation on x if a 4 is negative. 

Equating the coefficients of the cos 3/3z terms gives 

-0.89/1 



R = 



\l.l5A + 0.222 (l + ~) 1 



(110) 



Thus, with known R which is fully defined by the medium (a 2 and a 4 ) 
and the point of entry (xo) for the ray, one can compute A and 0. The 
plots of Figs. 28-30 show the interrelations between A, R, and /? as 
given by (108) and (110). For R « 1, A S -fl/4 and jS S Va 2 (1 + 
1.522/2). The principal effect of a small a 4 term in the index variation 
is to change the period of the ray oscillation. This is illustrated in Fig. 31 
for several values of R compared to the R = (square law) case. Posi- 
tive R means a 4 has a focusing effect, and the ray period is shortened. 
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Fig. 28— R vs A according to (110). 
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Fig. 29 — Normalized phase constant 0/v«2 vs A according to (108) and 
(110). 
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Fig. 30 — Normalized phase constant 0/y/a-i vs It according to (108) and 
(110). 
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Fig. 31 — Ray position path x/xo vs normalized distance for a particular 
non-square-law medium. 

With a renormalization of the abscissa, the curves of Fig. 31 are re- 
plotted in Fig. 32 to show that the cos 3/3z term is indeed small and the 
ray path differs little from a sine wave. 

However, the period 27r//3 does depend on the peak ray-path am- 
plitude (or ray-path slope at the axis), and hence the nice separation 
between the input-ray slope and input-ray position which was found in 
(10) does not exist for non-square law media. 

Because the ray period depends on the peak ray-path amplitude, a 
group of rays entering a non-square law medium at z = as in Fig. 33 
will fall out of step and at some large z one ray will be at a positive 
maximum when another is at a negative maximum. These rays can 
represent parts of a beam of light injected into the medium off-axis 
when the beam spot size is very large compared to a wavelength. This 
shows that the injected off-axis beam will spread out and occupy the 
region ±w t about the guide axis. However, the beam will never 
occupy any more than the region ±w t if it is injected parallel to the 
guide axis. 

One can obtain solutions analogous to the one given above if the 
input ray is on axis but at some slope x by noting 



x = Sc„ sin nftz 
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Fig. 32 — Ray position x/x vs ^ 2 /2tt for the medium of Fig. 31 showing basic 
similarity in shape of ray path. 




Fig. 33 — Ray position x vs distance z in a non-square-law medium illustrating 
the way an input beam breaks up as the wave propagates. 
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n = 1, 3, 5, etc. This follows from the same considerations which led to 
(47). Since we already have solved for the coefficients Ci and c 2 of (101), 
we can get the solution for (111 ) by transforming (101 ), letting x x (z ) = 
x (z — t/2/3) which yields 

%x = Ci sin 02 — c 2 sin 3/3z. (112) 

The parameters Ci and c 2 are related to x , the maximum of Xi , as in 
(105) and we define a new parameter B 

12 (*) 

B = ^ (113) 

which allows the ray slope to be written 

-r- = ci/3 cos &z — 3c2/3 cos 3/32 
dz 



= X\ { 1 + B sin (3z } cos fiz . 



- ./•;' {( 1 + V )cos|8z - ^cos3/fc) Ml '' 



The initial slope Xi is the ray slope at z = 0, and it is related to Ci and 
c 2 by 

*/ = cu3|l -H. (115) 

The previous interrelations between Ci , c 2 , A, /3, and i£ still hold. We 
seek a method for getting /3 and x , knowing only Zi , a^ and a 4 . To do 
so Xi, is rewritten in the form 

Xi = -TrVfid + A){1 + ft(1.5 + 0.3764 + 0.18754 2 )}*. (116) 

Hence, given x-l, (116) determines x/ s/aja* in terms of .ft (since R vs 
A is given in (110)); knowing ft, a 2 and a 4 , we can compute xo and 
using (109) and (108). Fig. 34 shows Xi s/ajai vs ft to facilitate this 
process. 

For an input ray with both slope and displacement, the proper match- 
ing to (101) with a suitable transformation can in principle be done, 
but has not been attempted. 

An approximate solution for the pure fourth-order medium 

n = n a {\ — \a& ) (117) 
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Fig. 34 — Normalized input ray slope Xi' vs R according to (116) and (110). 

in the equations which led to (108) and 



is obtained by letting a? 
(110). This yields 



/3 = x \/lA4a 4 



- = 23.3 



A = -0.165 
for the solution corresponding to (101) and 

/ 2 /— 

xi = x vo 4 



(118) 
(119) 
(120) 

(121) 



for the solution corresponding to (114). Equation (121) gives the ray 
slope at the axis x = 0. 

Better approximations would of course be obtained by including 
higher order terms in (101) but the large Ci/c 2 ratio produced by in- 
cluding only the cos 3/32 terms suggests that the cos 5/?2 term would 
have a negligible coefficient. 

We can combine (121) and (118) to express the ray period for the 
fourth-order medium in terms of the ray slope at the .r-axis; this gives 



= a 4 * Vl.44*i' • 



(122) 
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We can now compare the ray period from (122) with that previously 
derived, (85). We let x{ equal the characteristic ray angle denned by 
(81) with m = 0. Then the ray period from (122) becomes 

P Afl4 

The constant in (123) differs slightly from the value 12.3 found in (85) 
but the dependence on X and a\ is identical. 

IX. DISCUSSION AND ACKNOWLEDGEMENT 

The abstract contains a summing up review of the contents of this 
paper. One might add that a remarkably small step change in index is 
required to contain completely (for practical purposes) a light beam. 
As shown in Fig. 9, at X = 0.6328 n, a step change in index of refraction 
of a few parts in 10 6 is adequate for a beam radius of a = 0.419 mm, 
and this change need only be maintained from x = a to x ^ 2a where 
the energy is certain to be too small for that region to influence the 
wave propagation. 

The author would like to thank Mr. Tingye Li for the use of his 
computer program which was modified to make the computations repre- 
sented in Figs. 16-27. Without Mr. Li's previous work the author would 
not have included those figures which help justify the inferences leading 
to (59). Mrs. C. L. Beattie made the modifications and very effectively 
saw through all of the computations, for which the author is most 
appreciative. Mr. E. A. Marcatili on numerous occasions gave learned 
reactions to the newly developed ideas. 

APPENDIX 

We seek here to account for the expression (60) 

h -^W- (i24) 

As noted in connection with Figs. 11 and 12, in a medium with a 
continuous variation in index of refraction, the fields of the higher 
order modes extend farther from the axis than do the fields of the lowest 
order mode. The technique used here to determine the phase constant 
and characteristic ray angle for rather general media is to establish an 
equivalent width of medium in which the energy is completely confined 
— as in guides with perfectly conducting walls or with zero permittivity 
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walls. It is clear then that this equivalent width must be different for 
the lowest-order modes than for higher-order modes. 

When the author was easting about for a method of expressing this 
change, E. A. Mareatili pointed out that there is a characteristic radial 
distance at which the function describing the square-law medium's 
modes changes from an oscillating function to an exponentially decaying 
function. These functions are the parabolic cylinder functions and the 
value of x for the transition is* 



:, = 2 |/m + | 



(125) 



where m is the mode index. 

This was tried as an equivalent width, but was found to change much 
more rapidly at small m than the actual increase in extent of the field 
illustrated in Figs. 11 and 12. By examining the radii at which the field 
decreased to about 1 per cent of the peak value for various low-order 
modes, it was found that (124) represents the variation quite well. Thus, 
the general form (125), which is supported by the function theory for 
the square-law medium, was modified to fit actual known field width 
variations at small m. For large m, (124) and (125) do have the same 
variation. 

The method of defining the equivalent width, outlined in the body 
of the paper in connection with (59), causes media with higher than 
square-law variations in index to merge smoothly into the known be- 
havior of the step-change index variation including the gradual disap- 
pearance of the factor b. 
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